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Middle East respiratory syndrome coronavirus (MERS-CoV) is a new highly pathogenic human coronav- 
iruses that emerged in Jeddah and Saudi Arabia and has quickly spread to other countries in Middle East, 
Europe and North Africa since 2012. Up to 17 December 2014, it has infected at least 938 people with a 
fatality rate of about 36% globally. This has resulted in an urgent need to identify antiviral drugs that are 
active against MERS-CoV. The papain-like protease (PL pro ) of MERS-CoV represents an important antiviral 
target as it is not only essential for viral maturation, but also antagonizes interferon stimulation of the 
host via its deubiquitination activity. Here, we report the discovery that two SARS-CoV PL pro inhibitors, 
6-mercaptopurine (6MP) and 6-thioguanine (6TG), as well as the immunosuppressive drug mycophenolic 
acid, are able to inhibit MERS-CoV PL pro . Their inhibition mechanisms and mutually binding synergistic 
effect were also investigated. Our results identify for the first time three inhibitors targeting MERS- 
CoV PL pro and these can now be used as lead compounds for further antiviral drug development. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

Ten years after the outbreak of severe acute respiratory syn¬ 
drome (SARS) during 2002/2003 (Hilgenfeld and Peiris, 2013), a 
new highly pathogenic human coronavirus (CoV), Middle East 
respiratory syndrome coronavirus (MERS-CoV), has emerged in Jed¬ 
dah and Saudi Arabia; this virus has then quickly spread to other 
countries in the Middle East, Europe, and North Africa (Anderson 
and Baric, 2012; Chan et al., 2012; Zaki et al., 2012 ). The virus causes 
symptoms similar to SARS-CoV, but also leads to acute renal failure 
with the result that it is even more deadly than SARS (Eckerle et al., 
2013). As of 17 December 2014,938 people have been infected with 
MERS and this has led to 343 reported deaths (http://www.who.int/ 
csr/don/17-december-2014-mers/en/). SARS-CoV and MERS-CoV 
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belong to the genus, namely the Betacoronavirus (de Groot et al., 
2013). Bats are suspected to be their original reservoir as a few 
bat CoVs with high sequence similarity to SARS- and MERS-CoV 
have been identified (Drexler et al., 2014; Ithete et al., 2013). 
Recently, new evidence has suggested that dromedary camels are 
able to act as a zoonotic source of the virus (Haagmans et al., 
2014; Reusken et al., 2013). Limited human-to-human transmis¬ 
sion has also been evident in MERS outbreaks involving household 
and hospital contacts (Assiri et al., 2013). The significant number of 
clustered and sporadic cases that have originated from multiple 
sources indicate the potential for future outbreaks (de Groot 
et al., 2013; Ithete et al., 2013; Reusken et al., 2013). Such a situa¬ 
tion means that there is an urgent need to identify antiviral drugs 
against MERS-CoV in order to prepare for such future outbreaks. 

Proteases are one of the most prominent and useful drug targets 
in antiviral therapies; this includes proteases from human immu¬ 
nodeficiency virus (Wensing et al., 2010), hepatitis C virus 
(Hulskotte et al., 2012; Kwo et al., 2010), dengue virus (Steuer 
et al., 2011), and influenza (Zhirnov et al., 2011). Coronaviral pro¬ 
teases, specially the main protease (EC 3.4.22.69) and the papain¬ 
like protease (PL pro ) (EC 3.4.22.46), are considered to be suitable 
antiviral targets because they are responsible for the cleavage of 
nonstructural polyproteins (ppla and pplab) that are essential 
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for viral maturation. A large number of SARS-CoV main protease 
inhibitors have been identified and published (Anand et al., 

2003; Bacha et al., 2008; Verschueren et al., 2008; Wu et al., 
2006; Yang et al., 2003; Zhu et al., 2011). One type of protease 
inhibitors for SARS is the benzotriazole ester inhibitors; these tar¬ 
get the SARS-CoV main protease. These drugs have been identified 
as having inhibitory properties against the MERS-CoV main prote¬ 
ase (Kilianski et al., 2013). Unlike the drug’s broad-spectrum inhi¬ 
bition of the main protease inhibitors of SARS, initial screening of 
an existing SARS-CoV PL pro inhibitor, namely a benzodioxolane 
derivative BD-15 g, against MERS-CoV PL pro revealed no significant 
inhibition (Ghosh et al., 2010; Kilianski et al., 2013). The identifica¬ 
tion of this difference between MERS and SARS means that there is 
a need to further understand the structural and functional charac¬ 
teristics of MERS-CoV PL pro . 

In addition to proteolytic activity, which is similar to that of 
other coronaviruses, MERS-CoV PL pro has both deubiquitination 
(DUB) and ISGylation (ISG = interferon-stimulated gene) activity 
(Yang et al., 2014). The enzyme is able to deubiquitinate interferon 
regulatory factor 3, which prevents its nuclear translocation; this 
results in a suppression of interferon (3 production; the end result 
is immune suppression of host cells (Chen et al., 2007; Clementz 
et al., 2010; Yang et al., 2014; Zheng et al., 2008). Very recently, 
crystal structures of the MERS-CoV PL pro free enzyme and of 
MERS-CoV PL pro in complex with a ubiquitin (Ub) derivative have 
been published and these studies pinpoint the unique architecture 
of the oxyanion hole of MERS-CoV PL pro , which differs markedly 
from that of all other structurally characterized PL pro proteins 
(Bailey-Elkin et al., 2014; Lei et al., 2014). The unique features of 
its S3 and S5 subsites, as well as the flexible binding loop, provide 
potential drug targets within the active-site features that should 
allow structure-based drug design. Here, we report the discovery 
that two SARS-CoV PL pro inhibitors (Chou et al., 2008), 6-mercapto- 
purine (6MP) and 6-thioguanine (6TG), as well as the immunosup¬ 
pressive drug, mycophenolic acid, are able to inhibit MERS-CoV 
PL pro . Their inhibition mechanism and mutually binding synergistic 
effect were also investigated. These are the first three inhibitors of 
MERS-CoV PL pro to be described and our findings have implications 
with respect to future lead compounds and further drug 
development. 


2. Material and methods 

2.2. Expression and purification of MERS-CoV and SARS-CoV PL pro 

The protein induction and purification procedures for SARS- 
CoV, the E168R mutant of SARS-CoV PL pro and MERS-CoV PL pro 
have been described previously (Chou et al., 2012, 2014; Lin 
et al., 2014). The expression plasmid for MERS-CoV PL pro (GenBank 
accession number NC_019843.2; polyprotein residues 1484-1800) 
includes 6x His tag at the C-terminus (Lin et al., 2014). Briefly, the 
three expression vectors were separately transformed into Escheri¬ 
chia coli BL21 (DE3) cells (Novagen). These strains are incubated 
overnight at 20 °C and induced with 0.4 mM isopropyl-p-D-thioga- 
lactopyranoside. The cell pellets were resuspended in lysis buffer 
(20 mM Tris, pH 8.5, 250 mM NaCl, 5% glycerol, 0.2% Triton X- 
100, and 2 mM (3-mercaptoethanol), lysed by sonication and then 
centrifuged to remove the insoluble pellet. Next, the supernatant 
was incubated with 1-ml Ni-NTA beads at 4 °C for 1 h. After allow¬ 
ing the supernatant to flow through a column, the beads were 
washed with washing buffer (20 mM Tris, pH 8.5, 250 mM NaCl, 
8 mM imidazole, and 2 mM p-mercaptoethanol), and the protein 
was eluted with elution buffer (20 mM Tris, pH 8.5, 30 mM NaCl, 
150 mM imidazole, and 2 mM p-mercaptoethanol). The protein 


was then loaded onto a S-100 gel-filtration column (GE Healthcare) 
equilibrated with running buffer (20 mM Tris, pH 8.5, 100 mM 
NaCl, and 2 mM dithiothreitol). The purity of the fractions collected 
was analyzed by SDS-PAGE and the protein was concentrated to 
30 mg/ml using an Amicon Ultra-4 10-kDa centrifugal filter 
(Millipore). 

2.2. Deubiquitination (DUB) assay 

The DUB assay was carried out as described previously (Chou 
et al., 2008, 2014; Lin et al., 2014). Briefly, the fluorogenic substrate 
Ub-7-amino-4-trifluoro-methylcoumarin (Ub-AFC) (Boston Bio- 
chem) at 1 pM was incubated without or with the chemical com¬ 
pounds in 50 mM phosphate pH 6.5 for 3 min before the addition 
of the MERS-CoV or SARS-CoV PL pro of 0.17 pM. The enzymatic 
activity at 30 °C was determined by continuously monitoring, 
using fluorescence emission and excitation wavelengths of 350 
and 485 nm, respectively, in a PerkinElmer LS 50B luminescence 
spectrometer (USA). 

2.3. Steady-state kinetic analysis 

The fluorogenic peptidyl substrate, Dabcyl-FRLKGGAPIKGV- 
Edans, was used to measure the enzymatic activity of MERS-CoV 
and SARS-CoV PL pro , as well as the E168R mutant of SARS-CoV 
PL pro , as described previously (Chou et al., 2008; Lin et al., 2014). 
Specifically, the enhanced fluorescence emission upon substrate 
cleavage was monitored, using excitation and emission wave¬ 
lengths of 329 and 520 nm, respectively, in a PerkinElmer LS 50B 
luminescence spectrometer. Fluorescence intensity was converted 
to the amount of hydrolyzed substrate using a standard curve 
drawn from the fluorescence measurements of well-defined con¬ 
centrations of the Dabcyl-FRLKGG and APIKGV-Edans peptides in 
a 1:1 ratio. This approach also corrects for the inner filtering effect 
of the substrate. For the inhibition studies, the reaction mixture 
contained 4-50 pM of peptide substrate with 0-50 pM 6MP or 
6TG in 50 mM phosphate pH 6.5 or 4-50 pM of peptide substrate 
with 0-500 pM mycophenolic acid in 50 mM phosphate pH 6.5, 
all in a total volume of 1 mL. After the addition of the enzyme to 
the reaction mixture, the increase in fluorescence was continu¬ 
ously monitored at 30 °C. The increase in fluorescence was linear 
for at least 3 min, and thus the slope of the line represented the ini¬ 
tial reaction velocity ( v ). 

The inhibition data for 6MP or 6TG were found to best fit a com¬ 
petitive inhibition pattern according to Eq. (1): 

v = k cat [E][S]/((l + [I\/K[ S )K m + [S]) (1) 

while inhibition data for mycophenolic acid were found to best 
fit a noncompetitive inhibition pattern according to Eq. (2): 

v = k C3t {E\[S]/((\ + {I]/K is )(K m + {S\)) (2) 

in which /< cat is the rate constant, [E], [S] and [I] denote the 
enzyme, substrate and inhibitor concentrations, and K m is the 
Michaelis-Menten constant for the interaction between the pep¬ 
tide substrate and the enzyme. K is is the slope inhibition constant 
for the enzyme-inhibitor complex. The program SigmaPlot 12 
(Systat Software Inc., USA) was used for the data analysis. 

2.4. Multiple inhibition assay 

To characterize the mutual effect of the inhibitors, the activity 
of PL pro was first measured with or without either 6MP (0 and 
8 pM) or 6TG (0 and 5 pM) in the presence of various concentra¬ 
tions of mycophenolic acid (0-180 pM). Then the activity of PL pro 
was measured with or without 6MP (0 and 8 pM) or mycophenolic 
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acid (0 and 150 jiM) at various NEM concentration (0-50 jiM). The 
peptidyl substrate and the enzyme concentrations were held at 15 
and 1 pM, respectively. Data obtained from the multiple inhibition 
reactions were fitted to the Eq. (3): 

v = Vo/(l + [I\/K, + [J]/Kj + [I][J\/«K,Kj) (3) 

where v is the initial velocity in the presence of both inhibitors, [I] 
and [/] are the concentrations of the two inhibitors, v 0 is the velocity 
in the absence of inhibitors, and Kj are the apparent dissociation 
constants for the two inhibitors, and a is a measure of the degree of 
interaction of the two inhibitors (Yonetani and Theorell, 1964). 

2.5. Computer modeling of PL pro in a complex with 6MP or 
mycophenolic acid 

The crystal structure of MERS-CoV PL pro (PDB code: 4PT5) was 
used as the template. As described previously (Chen et al., 2009; 
Chou et al., 2008), docking was performed using DS Modeling 1.7 
software (Accelrys). Before docking, the space near the catalytic 
triad (Cyslll-His278-Asp293) was chosen for the docking. The 
structures of 6MP and mycophenolic acid were then created and 
separately specified as input ligands. During the docking, chemical 
flexibility of ligands was allowed, while the grid extension from 
the site and the nonbonded cutoff distance were set as 3.0 and 
10.0 A, respectively. The docking results were calibrated using 
Monte Carlo trials, in which possible poses of ten were compared. 
After energy minimization, the pose with the lowest binding 
energy was chosen for the binding model. 

3. Results 

3.1. Identification of inhibitors of MERS-CoV PL pro 

Previous studies have indicated that two thiopurine analogs, 
6MP and 6TG, are able to inhibit SARS-CoV PL pro and that their 
binding sites are close to the catalytic triad (Chou et al., 2008, 
2014). Recently, Hilgenfeld and colleagues determined the struc¬ 
ture of MERS-CoV PL pro by X-ray crystallography (Lei et al., 2014) 
and confirmed that the SARS and MERS PL pro proteins have similar 
architecture; albeit they only show 30% sequence identity and 50% 
similarity. These findings led us to speculate that thiopurine ana¬ 
logs might be able to inhibit the MERS-CoV PL pro because the 
two PL pro possess a similar catalytic triad. In addition, another 
pre-existing drug, mycophenolic acid, has been identified as an 
effective anti-MERS-CoV agent (Chan et al., 2013), although details 
of the drug’s antiviral effect mechanism remain unclear. In the 
present study, we have evaluated the inhibitory effect of these 
three compounds on MERS-CoV PL pro . Using peptide substrates 
and Ub-AFC, we found that 6MP, 6TG, and mycophenolic acid are 
all independently able to inhibit effectively the proteolytic activity 
and deubiquitination of MERS-CoV PL pro (Fig. 1 and Supplementary 
Fig. 1). However, mycophenolic acid was found to be unable to 
inhibit SARS-CoV PL pro even at a concentration of 1 mM; neverthe¬ 
less it shows a dose-dependent inhibitory effect on the E168R 
mutant of SARS-CoV PL pro (Fig. 1). As the residue 168 in MERS- 
CoV PL pro is also an arginine (Chou et al., 2014), this suggests that 
mycophenolic acid is able to recognize the positive charge of resi¬ 
due 168 of PL pro . Among the three inhibitors, 6MP and 6TG have 
IC 50 values of 12.4-26.9 pM, while mycophenolic acid has an IC 50 
of 222.5-247.6 pM ( ^able 1). As a comparison, NEM, which is a 
commonly used cysteine protease inhibitor that acts by covalently 
modifying the active-site Cys of proteases (Chou et al., 2008), was 
also verified to be an inhibitor of MERS-CoV PL pro , this time with an 
IC 50 of 45 pM (Table 1). Compared with NEM, 6MP and 6TG are 
more effective inhibitors against the MERS-CoV PL pro , while 


mycophenolic acid is a less effective inhibitor against the 
MERS-CoV PL pro . 

3.2. Inhibition mechanism 

To understand the kinetic mechanisms involved in the interac¬ 
tion of these compounds with MERS-CoV PL pro , the proteolytic 
activity of the enzyme was measured using a series of peptidyl 
substrate concentrations and at various inhibitor concentrations. 
The inhibition data were then globally fitted to all possible kinetic 
models (competitive, noncompetitive and uncompetitive). A sim¬ 
ple competitive inhibition pattern best described the experimental 
data obtained for 6MP (Fig. 2A) and for 6TG (Fig. 2B) with the dou¬ 
ble reciprocal plot showing all lines intercepting they-axis. In con¬ 
trast, a noncompetitive inhibition pattern best described the 
inhibition caused by mycophenolic acid with the double reciprocal 
plot showing all lines intercepting the x-axis (Fig. 2C). The best fit 
of Eq. (1) to the experimental data obtained using 6MP and 6TG 
and of Eq. (2) to the experimental data obtained using mycophen¬ 
olic acid yielded inhibition constant (7C is ) values of 14.3, 9.1 and 
263.7 pM, respectively (Table 2). As summarized in Table 2, the 
Michaelis constant (J( m ) for peptidyl substrate hydrolysis in the 
absence or presence of inhibitors ranged from 14.1 to 21.9 pM 
and the k cat values ranged from 0.23 to 0.26 min -1 . The I< m for 
the peptidyl substrate was 1.5-fold higher than of the I< is for 6MP 
and 6TG, while the I< m for the peptidyl substrate was 17.3-fold 
lower than of the I< [s for mycophenolic acid. 

3.3. Binding synergy analysis 

As a result of different inhibition mechanisms at work when the 
thiopurine analogs and mycophenolic acid are present, this raised 
our interest in the characterization of the mutual binding effects of 
the two kinds of compounds; as a result we carried out multiple 
inhibition assays (Chou et al., 2014; Yonetani and Theorell, 
1964), which are shown in Fig. 3. The experiments were carried 
out by measuring the initial reaction velocities of PL pro at various 
concentrations of mycophenolic acid with either 6MP or 6TG pres¬ 
ent. The substrate concentrations were held constant at subsatura¬ 
tion levels, which allow both inhibitors to bind to the enzyme 
(Benson et al., 2008). By fitting the data to Eq. (3), both results 
showed that the lines intersect above the x axis using the Yonetan- 
i-Theorell plot and that the a value for mycophenolic acid in com¬ 
bination with 6MP is 0.73 and for mycophenolic acid in 
combination with 6TG is 0.19 (Fig. 3A and B). These findings show 
that the binding of mycophenolic acid and each thiopurine 
together result in a synergistically enhancement of the inhibition 
effect of the pairs of drugs on MERS-CoV PL pro . On the other hand, 
when the same analysis was carried out using various concentra¬ 
tions of NEM with either 6MP (Fig. 3C) or mycophenolic acid 
(Fig. 3D), the results were different. Interestingly, the data points 
at various NEM and 6MP formed parallel lines (Fig. 3C) and the 
best-fit oc value is one that is nearly infinite (2.7 x 10 7 ). These find¬ 
ings indicate that NEM and 6MP bind in a mutually exclusive fash¬ 
ion (Copeland, 2000). Based on these results and that fact that NEM 
covalently binds to the cysteine residue in the active site, this sup¬ 
ports the hypothesis that 6MP can also bind to the enzyme’s cata¬ 
lytic triad (Cyslll-His278-Asp293) and that this results in its 
inhibitory effect. On the other hand, the binding synergy assay 
using various concentration of NEM in the presence of mycophen¬ 
olic acid showed that the best-fit lines intersected above the x axis 
in the Yonetani-Theorell plot and had an a value of 0.22 (Fig. 3D). 
These findings indicate that the binding of NEM and mycophenolic 
acid are similar to the findings for 6MP/6TG and mycophenolic 
acid, namely that NEM and mycophenolic acid are also able to 
synergistically increase MERS-CoV PL pro inhibition. 
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Fig. 1. Inhibition of coronaviral PL pro proteolytic activity by thiopurines, NEM, or mycophenolic acid. The proteolytic activity of MERS-CoV and SARS-CoV PL pro , as well as the 
E168R mutant of SARS-CoV PL pro in the presence of 6MP (8-50 pM), 6TG (6-50 pM), NEM (50-500 pM) or mycophenolic acid (50-2000 pM) were measured, respectively. The 
concentrations of peptidyl substrate (Dabcyl-FRLKGGAPIKGV-Edans) was 20 pM, while the concentration of MERS-CoV PL pro , SARS-CoV PL pro , and its E168R mutant was 1, 
0.1, and 0.3 pM, respectively. MPA: mycophenolic acid. 


Table 1 

Structure and IC 50 of compounds against coronaviral PL pro . 


Name 


Structure 


6-Mercaptopurine (6MP) 


6-Thioguanine (6TG) 


N-ethylmaleimide (NEM) 


Mycophenolic acid 



IC 50 (pM) a 

Peptide cleavage DUB activity 

MERS-CoV PL pro 

26.9 ±7.5 25.8 ±5.0 


24.4 ± 4.3 


12.4 ± 1.9 


45.0 ± 10.1 


ND b 


MERS-CoV PL pro 

247.6 ±55.8 222.5 ± 57.8 

SARS-CoV PL pro 

No inhibition No inhibition 

SARS-CoV PL pro E168R 

260.1 ±21.0 ND b 


a The peptide substrate (Dabcyl-FRLKGGAPIKGV-Edans) at 20 pM or Ub-AFC at 0.5 pM were used for the IC 50 assay, respectively. All the assays were repeated several 
times to ensure reproducibility. 
b Not determined. 


3.4. Putative binding site analysis 

To gain further insights into the inhibition mechanism, we per¬ 
formed docking experiments using 6MP and mycophenolic acid 
and MERS-CoV PL pro . Both compounds are able to inhibit the 
deubiquitination and proteolytic activity of the enzyme, suggesting 
that their binding may be close to the binding site for the Ub tail 
(Leu73-Arg74-Gly75-Gly76), not that of the Ub core (residues 
1-72) (Chou et al., 2014). Therefore, the space near the catalytic 
triad (Cyslll-His278-Asp293), from the putative S4-S2' subsites, 
was chosen as the docking site. The best docking scores for 6MP 
and mycophenolic acid binding to MERS-CoV PL pro were 15.1 and 
31.8, respectively. Moreover, the binding energy of 6MP and myco¬ 
phenolic acid to PL pro were -18.4 and -34.6 kcal/mol, respectively. 
Interestingly and importantly, the docked sites for the two 
compounds are quite different (Fig. 4). 6MP fits well into the 
cavity near the catalytic triad of PL pro and it is able to make 


hydrogen-bonding interactions with residues Asnl09, Gly277 and 
His278. This putative binding site is close to that found for SARS- 
CoV PL pro (Fig. 4, carbon atoms colored in gray) (Chou et al., 
2008). After drug binding, these hydrogen-bonding interactions 
with the enzyme may interfere with proton transfer from the thio 
group of Cyslll to the imidazole ring of His278; this is likely to 
bring about a failure of catalysis. In contrast, mycophenolic acid 
was found to show hydrogen-bonding interactions with residues 
Aspl64 and Argl68 and to have hydrophobic contact with residues 
Prol63, Val276 and Phe269 (Fig. 4). In terms of comparison with 
SARS-CoV PL pro , this pocket is close to the residues Glyl64 and 
Glul68 of SARS-CoV PL pro . However, when SARS-CoV PL pro was 
used for the docking of mycophenolic acid, this failed. This finding 
is able to explain why mycophenolic acid shows an inhibitory 
effect on MERS-CoV PL pro , but not on SARS-CoV’s. Furthermore, 
when we compared our findings with the structures of coronaviral 
PL pro -Ub complex (Bailey-Elkin et al., 2014; Chou et al., 2014), the 
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Fig. 2. Inhibition of MERS-CoV PL pro by 6MP (A), 6TG (B), or mycophenolic acid (C). 
The proteolytic activity of PL pro was measured in the presence of different peptide 
substrate concentration (4-50 pM) and various inhibitor concentrations (0-50 pM 
for 6MP or 6TG and 0-500 pM for mycophenolic acid). The solid lines in panel A and 
B are the best fit by the competitive inhibition (Eq. (1)) and those in panel C 
are the best-fit by the noncompetitive inhibition (Eq. (2)). The R sqr are 0.978, 0.979, 
and 0.988, respectively. All assays were repeated several times to ensure 
reproducibility. The kinetic parameters such as I< m , k cat and K is from the best fit 
were shown in Table 2. 


location of the mycophenolic acid docking site was found to be 
close to the putative S3-S4 subsites of PL pro . This suggests that 
mycophenolic acid is able to interfere with the binding or orienta¬ 
tion of the substrate P3-P4 residues. Overall, the fact that we found 
distinct binding sites for the compounds studied here conforms 
with their synergistic inhibition of MERS-CoV PL pro . 

4. Discussion 

6MP and 6TG have been used clinically to treat cancer since the 
1950’s. Up to the present, hundreds of 6MP and 6TG related studies 
have been published every year, including metabolism, pharmaco¬ 
kinetics and dosage as well as their use in various therapies that 
target cancers, rheumatoid arthritis, systemic lupus erythematosus 
and Crohn’s disease (Karran and Attard, 2008). 6TG has recently 
been proposed as a means of treating a wide variety of cancers that 
have a high frequency of homozygous deletion of the gene encod¬ 
ing methylthioadenosine phosphorylase, which include leukemia, 
glioblastoma, lung cancer, pancreatic, squamous cell carcinoma 
and others (Munshi et al., 2014). In addition, 6MP and its pro-drug, 
azathioprine, have also been used for many years as an immuno¬ 
suppressant when treating organ transplant recipients (Karran 
and Attard, 2008). Furthermore, our previous studies have sug¬ 
gested that 6MP and 6TG are reasonably selective inhibitors of 
PL pro and do not inhibit other cysteine proteases such as the coro- 
naviral main protease, cathepsins B, K, L, and S, and papain (Chou 
et al., 2008). Very recently, 6MP has been identified to have a pro¬ 
tective effect against influenza A/WSN/1933 (H1N1) virus, albeit 
with a high EC 50 of 26.5 pM and a low selectivity index of 4 
(Chan et al., 2013). These studies suggest that thiopurines have a 
broad-spectrum of activity and are likely to be a suitable progeni¬ 
tor for the development of anticancer and antiviral drugs. 

An immunosuppressant, mycophenolic acid, has been shown to 
have broad activity in vitro and/or in animal models against differ¬ 
ent viruses, including Chikungunya virus (Khan et al., 2011), den¬ 
gue virus (Diamond et al., 2002), Japanese encephalitis virus 
(Sebastian et al., 2011), West Nile (Morrey et al., 2002), yellow 
fever virus (Leyssen et al., 2005), and possibly hepatitis B virus 
(Pan et al., 2012). Recently, Yuen and colleagues in Hong Kong have 
found that a combination of mycophenolic acid and interferon pib 
is able to synergistically inhibit the replication of MERS-CoV in 
Vero cells; albeit their pharmacological mechanism of action 
remain unknown (Chan et al., 2013). Here we investigated whether 
mycophenolic acid and 6MP or 6TG are able to synergistically inhi¬ 
bit MERS-CoV PL pro . The present studies provide the molecular 
basis for the design of more effective PL pro inhibitors against MERS. 

GRL-0617, a highly effective SARS-CoV PL pro inhibitor whose 
binding site is in the S3-S4 site, also has a noncompetitive inhibi¬ 
tion mechanism (Lee et al., 2013). As there are significant differ¬ 
ences in the S3 subsite, as well as in the BL2 loop, between the 
SARS-CoV and MERS-CoV PL pro s, it is unlikely that GRL0617 and 
its derivatives will also inhibit MERS-CoV PL pro (Lei et al., 2014). 
However, the discovery of the fact that different compounds are 


Table 2 

Kinetic parameters of 6MP, 6TG, and mycophenolic acid inhibition of MERS-CoV PL pr0 . 


PL pro in 

Krn (nM) 

kcat (min ^ 

k C at/K m (10 3 min 1 pM : ) 

K is (|.iM) 

pH 6.5 a 

19.2 ±2.6 

0.24 ± 0.01 

12.5 ± 1.9 


6MP b 

21.9 ±3.3 

0.28 ± 0.02 

12.8 ±2.0 

14.3 ±1.6 

6TG b 

14.1 ±2.3 

0.23 ± 0.02 

16.3 ±3.0 

9.1 ± 1.2 

Mycophenolic acid 

15.2 ± 1.1 

0.26 ± 0.01 

17.1 ± 1.4 

263.7 ± 13.8 


a The values were from our previous studies (Lin et al., 2014). 

b In the presence of 6MP or 6TG, the steady-state kinetic data were determined by the competitive inhibition model (Eq. (1)) (Fig. 2A and B). 
c In the presence of mycophenolic acid, the kinetic data were determined by the noncompetitive inhibition model (Eq. (2)) (Fig. 2C). 
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[Mycophenolic acid] (uM) 


[NEM] (uM) 


Fig. 3. Mutual effect of the coronaviral PL pro inhibitors. Two inhibitor patterns for mycophenolic acid with either 6MP (A) or 6TG (B) present, and those for NEM with either 
6MP (C) or mycophenolic acid (D) present were performed. The concentrations of peptidyl substrate and enzyme were 15 and 1 pM, respectively. The points are the reciprocal 
of the experimental velocities, and the lines are the best fit of the data to Eq. (3). The results suggest that ot value of the four experiments (A-D) is 0.73, 0.19, 2.7 x 10 7 and 
0.22, respectively. 
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Fig. 4. Docking of 6MP and mycophenolic acid with MERS-CoV PL pro . The docking 
for 6MP (carbons colored by orange) and mycophenolic acid (magenta) with MERS- 
CoV PL pro (cyan; PDB code: 4PT5) were performed with DS modeling 1.7 software 
(Accelrys). The docked models and that for 6MP with SARS-CoV PL pro (gray) are then 
overlaid for comparison. The dotted lines show the hydrophilic interactions 
between PL pro and inhibitors. 6MP is able to be targeted to the catalytic triad of 
MERS-CoV PL pro . Instead, mycophenolic acid is located far from the catalytic triad 
and shows hydrogen bonding interactions with the residues Aspl64 and Argl68, 
which are Glyl64 and Glul68 in SARS-CoV PL pro . According to the structure of 
SARS-CoV PL pro in complex with Ub, the docked site of mycophenolic acid should be 
near the putative S3-S4 subsites, whereas 6MP may occupy SI’ subsite (Chou et al., 
2014). The figure was produced using PyMol (http://www.pymol.org). (For inter¬ 
pretation of the references to color in this figure legend, the reader is referred to the 
web version of this article.) 


able to bind to a similar region within MERS-CoV PL pro is to be 
expected. Based on the selective inhibitory effect of mycophenolic 
acid on MERS-CoV PL pro and on SARS-CoV PL pro E168R mutant 
(Fig. 1), we can conclude that mycophenolic acid is able to recog¬ 
nize the positive charge of residue 168, which is near the S3-S4 
subsites (Bailey-Elkin et al., 2014; Chou et al., 2014). Indeed, our 
docking experiments suggest that the binding site of mycophenolic 
acid is quite close to the S3-S4 subsites. Moreover, because myco¬ 
phenolic acid cannot inhibit SARS-CoV PL pro , this finding may be 
used to explain why mycophenolic acid is ineffective against 
SARS-CoV using an animal model (Barnard et al., 2006). Although 
the affinity of mycophenolic acid with respect to MERS-CoV PL pro 
is weak, this drug can be used as a lead compound, which will 
allow further optimization to enhance potency. On the other hand, 
another study has shown that mycophenolic acid is a very effective 
inhibition of MERS-CoV replication in Vero cells, with an EC 50 of 
0.27 pM (Chan et al., 2013). These conflicting findings at such a 
low dose suggest that mycophenolic acid may have one or more 
other inhibitory targets or that there are additional effects against 
MERS-CoV. 

There have been up to the present no potent inhibitors that tar¬ 
get MERS-CoV PL pro . In the present study we found one noncom¬ 
petitive inhibitor (mycophenolic acid) and two competitive 
inhibitors (6MP and 6TG) that target the MERS-CoV PL pro synergis- 
tically. As both thiopurines and mycophenolic acid are still used 
extensively in a range of clinical setting, particularly when decid¬ 
ing leukemia treatment and during immunosuppression, our study 
has emphasized the need to study in depth the in vivo activity of 
these drugs in combination. Further evaluation of the potential 
therapeutic effects of these commercially available drugs in terms 
of in vitro activity may be needed because this virus is able to 
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spread globally and continues to pose a significant threat to world 
health and the world’s economy. 
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